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Получено интегральное уравнение трехмерной контакт- 
ной задачи для ортотропного полупространства (9 неза- 
висимых упругих параметров в законе Гука), ядро кото- 
рого не содержит квадратур и зависит от решения ха- 
рактеристического бикубического уравнения. Рассмот- 
рено взаимодействие двух одинаковых симметрично 
внедряемых жестких штампов, имеющих форму эллип- 
тических параболоидов. При неизвестной области кон- 
такта для решения этой задачи использован метод не- 
линейных граничных интегральных уравнений Галано- 
ва, позволяющий одновременно определить область 
контакта и давления в этой области. Для отладки ком- 
пьютерной программы использовано точное решение 
для одного эллиптического штампа. При заданной осад- 
ке, форме основания и взаимной удаленности штампов 
для разных ортотропных материалов рассчитаны кон- 
тактные давления, области контакта и вдавливающие 
силы. Модель ортотропного тела применяется для опи- 
сания многих востребованных в технике и промышлен- 
ности материалов: сера, сегнетовая соль, вольфрамит, 
барит, древесина различных пород. 


Ключевые слова: теория упругости, контактные зада- 
чи, ортотропное полупространство, взаимодействие 
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Введение. Уравнения упругого равновесия и закон Гука для ортотропного тела описаны в монографии [1]. 
Примеры ортотропных материалов даны в работах [2, 3]. Интегральное уравнение (ИУ) трехмерной контактной зада- 
чи для ортотропного полупространства, ядро которого выражено через двукратный интеграл, и его точное решение 
для кругового штампа впервые было получено А. О. Ватульяном [4]. В работах [5, 6] предложен метод освобождения 
от квадратур в ядре ИУ для трансверсально изотропного полупространства, основанный на теории обобщенных функ- 
ций и применимый также для ортотропного полупространства. В результате существенно упрощается расчет и регу- 
ляризация ядра ИУ, что и позволяет применить для решения контактных задач метод Галанова [7]. Исследовались 
точные решения контактных задач [8, 9] и взаимодействие штампов [10] для трансверсально изотропного полупро- 
странства. Цель настоящего исследования — изучить взаимодействие двух одинаковых штампов на ортотропном по- 
лупространстве. 





* Работа выполнена по гранту РФФИ 15-01-00331. 
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Контактная задача. В декартовых координатах рассмотрим ортотропное упругое полупространство 520. Оси 
упругой симметрии совпадают с осями координат. Закон Гука в прямой форме (выражения напряжений через дефор- 
мации) включает 9 независимых упругих параметров с; (#=1, 2, 3, 4, 5, 6), с1›, саз, Сэз [1]. Пусть при <=0 в полупро- 
странство внедряются два одинаковых абсолютно жестких штампа (эллиптические параболоиды, вершины которых 
расположены на оси х), основания которых описываются функциями 


2 2 
Р-Я, ЩЬЮ, 
2 2» 

Пусть задача симметрична относительно оси у. Штампы вдавливаются без перекоса одинаковыми силами Р, 
испытывая осадку 6. При заданных упругих параметрах, величинах №1, К», Й и осадке 6 требуется определить области 
контакта 62., контактное давление 4(х,у)=с.(х,у,0)/сзз и силу Р. 

С учетом симметрии задачи 4(-х.у)=а(х,у). Тогда ИУ задачи можно свести к ИУ на одном участке. После замен 
х, =х-П, 4, (х,, У) = 4(%, 9), ©, О, 


это ИУ можно переписать в форме (звездочки далее опускаем): 





[[челк(л- у. 5=®+ К - ,5+х+ 28) мым = 2448-70 У), 55), (1) 
о 
2 2 
т в >К.. 


Ядро ИУ (1) представимо в форме свободной от квадратур: 





а0  шиоиз (и + 2) (М5 + из)(из + Ш) _ Е(51,52) @) 
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2,2 : 
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ру = ао [Аут + 2(2у6 +Уз —9778)51 52 +4р53 |, 


К(151, 752) = 


ро = -—(у55ё +452 ув (151 + 252) + (А-2узув +2у6777з)51 52]. 


Здесь ша, и2, из — корни уравнения 
арий +а›и“ +ади? +ав =0, Вен, >0, в =у4у5, (3) 
а› =-[у4(А, — 2577) -+576 51 — [75 (А> —29478)+7476152, 
ад =[у6 (А, —29577)+117415154 + [У (А› 29478) +757415]52 + 





А 
+251 55 [5+9 (44 +185 +17) +1415 узтув — 1748 = 721577 — 1316], 


ав =-(у55? +1452 [ув (191 +1252) (улт2 — 29376 -73)5152], 
— . -а ый = 
У; ат 7=1,2,4,5,6, уз = с11с33, У7 = сс, Уз = С23Сзз, То =7Уз + Тб, 


А, — 11 —у7, Д> =7> —уз, А=сз Че ст , пт =1,2,3. 








При вычислении ядра (2) в каждой точке приходится решать новое кубическое характеристическое уравнение 
(по формулам Кардано), получающееся из уравнения (3). 

В табл. 1 даны значения безразмерных параметров 7; (=1, 2, ..., 8) (3) для ряда материалов [2, 3]. 

Для решения ИУ (3) при условии 4(х,у)=0, (х,у)ЕОО, используем метод нелинейных граничных ИУ типа 
Гаммерштейна, позволяющий одновременно определить область контакта и контактное давление. Суть метода 
подробно изложена в работах [7, 10]. 

Предположим, что область контакта © в ИУ (3) целиком содержится в прямоугольнике 


= ха, УЮЫ}. 
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Таблица 1 
Значения характеристик (безразмерные) 
Материал у! 2 уз 4 5 Ув у7 Уз 

Топаз, АБ(Е, ОН). 0,956 1,183 | 0,427 | 0,366 0,451 | 0,444 | 0,288 | 0.298 
Ангидрит, Саб О. 0,838 1,652 | 0,147 | 0,290 0,237 | 0,0827 | 0,136 | 0,283 
Сера, 5 0,497 0,424 | 0,275 | 0,0890 | 0,180 | 0,157 | 0,354 | 0,329 

Барит, ВабО. 0,832 0,757 | 0,428 | 0,112 0,251 | 0,240 | 0,283 | 0,266 
Целестин, 57504 0,812 0,825 | 0,601 | 0,105 0,217 | 0,207 | 0,470 | 0,481 
ны Вы 0,758 0,719 | 0,358 | 0,216 0,271 | 0,0871 | 0,341 | 0,295 

Сегнетова соль 

о 0,687 1,027 | 0,380 | 0,361 | 0,0865 | 0,264 | 0,313 | 0,394 
Ясень белый 0,162 0,104 | 0,0628 | 0,0702 | 0,0994 0 0,0878 | 0,0736 
Береза желтая 0,109 | 0,0698 | 0,0484 | 0,0637 | 0,0693 | 0,0159 | 0,0686 | 0,0532 
Дуб красный 0,181 | 0,0958 | 0,0559 | 0,0760 | 0,0835 0 0,0885 | 0.0631 
Грецкий орех черный 0,142 | 0,0749 | 0,0575 | 0,0558 | 0,0765 | 0,0189 | 0,107 | 0,0759 





Лиственница западная 0,0915 0,0753 0,0322 0,0690 0,0630 0,0070 0,0414 0,0312 





Сосна широкохвойная 0,117 0,0633 0,0402, 0,0600 0,0710 0,0120 0,0536 0,0312 



































Прямоугольник 5 покроем равномерной сеткой из т узлов с шагами й! по оси х и й> по оси у. При расчете 
значений ядра в этих узлах его особенности сглаживались по формулам 
2 2 2 и т 4 
(Е) ШУ 6-0 +6 Е. (4) 
Регуляризация (4) обеспечивает сходимость метода и отладку программы, давая хорошее совпадение с точным 
решением для одного эллиптического параболоида, которое имеет вид 


52 у? ый у? 
а(х, у) = ао 2 О: м =, 























а? 
2 2 
х 2 
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а’ Ь э 
о 
где 
5= а Бао | с(ф)аф : с(ф) = Е(соз? ф,зш? ф), 
и \/а? 05? ф+Ь? 5? Ф 
2 2 
К. _ С , а 4 Ь = 5, (6) 
в а 2 28 
Е. 2 
с= | с(ф)соз? фф = | Ри вы 
[2 05? ф+зт? Фр? ь [52 с05? ф+зш? ФР? ь Ь г 


0 0 

При заданных величинах 6, №1, К› сперва из формулы (6) для К/К! следует найти 5, а затем определить полуоси 
эллипса контакта а и В из следующей формулы (6). Величина ди далее находится из первой формулы (6). 

Отметим, что даже при вдавливании кругового штампа (К =) область контакта в решении (5) будет в общем 
случае эллиптической. 

Численный анализ. При расчетах брали @0=Бо (прямоугольник 5 — квадрат). Число узлов бралось от 81 до 289 
(при этом наибольшая погрешность вычислений наблюдается вблизи границы области контакта). 

Введем безразмерные обозначения (штрихи далее опускаем; контактное давление в (3) уже является 
безразмерным) 
2, у 5 О Юр, и ф. — 
а а а 0 Чо 0 0 0 
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В табл. 2 приведены значения 4и=4(0,0) и Р для случая вдавливания круговых параболоидов (А!=К›=0,5) при 
6=1. Случай й=0 соответствует вдавливанию только одного штампа. 

Как показывают расчеты, давление и сила возрастают с ростом осадки штампов. При сближении штампов (с 
уменьшением значения й) проявляется сложный характер взаимодействия штампов: функция распределения 
контактных давлений становится несимметричной, уменьшаются значения давления и силы. 


Таблица 2 
Значения давления и силы 




































































= й=5 п=,1 
Материал 
Чо Р Чо Р Чо Р 

Топаз 0,886 0,922 0,873 0,881 0,836 0,763 
Ангидрит 0,758 0,788 0,748 0,757 0,720 0,661 
Сера 0,514 0,531 0,506 0,504 0,481 0,433 
Барит 0,636 0,656 0,626 0,623 0,595 0,534 
Целестин 0,557 0,575 0,548 0,546 0,521 0,468 
Вольфрамит 0,695 0,724 0,685 0,691 0,656 0,599 
Сегнетова соль 0,652 0,663 0,645 0,640 0,622 0,572 
Ясень 0,433 0,450 0,426 0,429 0,408 0,371 
Береза 0,393 0,409 0,387 0,391 0,371 0,339 
Дуб 0,425 0,442 0,418 0,422 0,401 0,366 
Орех 0,392 0,407 0,386 0,388 0,369 0,335 
Лиственница 0,388 0,404 0,383 0,387 0,367 0,337 
Сосна 0,391 0,406 0,385 0,388 0,368 0,336 








Материалы из табл. 1 могут быть разделены на 2 типа. К первому типу относятся материалы, для которых вы- 
полнено следующее: если к поверхности полупространства из этого материала в начале координат приложена сосре- 
доточенная сила (4(х.у)=6(х)6(у)), то нормальное перемещение поверхности в точке х=0, у=1 будет больше, чем в точке 
х=1, у=0 (ангидрит, сегнетова соль и лиственница). Для остальных материалов (второго типа) больше будет переме- 
щение в точке х=1, у=0. В случае й=о (при внедрении одного кругового параболоида) область контакта вытягивается 
вдоль той оси, на которой меньше нормальные перемещения точек равноудаленных от начала координат, где прило- 
жена нормальная сосредоточенная сила. 

Как показывают расчеты области контакта, сделанные для сегнетовой соли и барита (К=К›=0,5, 6=1), в срав- 
нении со случаем одного штампа (й=оо) взаимодействие близко расположенных штампов (й=1,1) приводит к уменьше- 
нию площади области контакта под каждым штампом (при заданной осадке). 

Выводы. Представление ядра ИУ в форме, свободной от квадратур, позволяет эффективно исследовать взаи- 
модействие штампов на ортотропном полупространстве. Аналогично может быть исследован случай, когда вершины 
штампа расположены на оси у. Развитый метод также может быть применен для случая конических, пирамидальных и 
других штампов. С практической точки зрения наиболее важными представляются расчеты на контактную прочность 
для многочисленных пород древесины, которые обладают ортотропной структурой. 
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